Fourth generation quark mass limits in CKM-element space 
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We present a reanalysis of CDF data to extend limits on individual fourth-generation quark masses 
from particular flavor-mixing rates to the entire space of possible mixing values. Measurements from 
CDF have set individual limits on masses, m' b and m' t , at the level of 335-385 GeV assuming specific 
and favorable flavor-mixing rates. We consider the space of possible values for the mixing rates and 
find that the CDF data imply limits of 290 GeV and greater over a wide range of mixing scenarios. 
We also analyze the limits from the perspective of a four-generation CKM matrix. We find that 
present experimental constraints on CKM elements do not suggest further constraints on fourth- 
generation quark masses. 
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I. INTRODUCTION 

A simple modification of the standard model is the ad- 
dition of a fourth sequential generation of fermion dou- 
blets. This natural extension of the standard model [1] 
may trigger dynamical electroweak symmetry break- 
ing [2] without a Higgs boson, and so address the hierar- 
chy problem. The new heavy fermions may have Yukawa 
couplings so large that they become strong. The subse- 
quent strong dynamics may lead to a composite of fourth 
generation fermions performing the role of the Higgs [3- 
5]. 

Recent searches by the CDF Collaboration for direct 
production of the fourth generation quarks, denoted t' 
and b' for the up- and down-type, found m t / > 335 
GeV[6] and my > 385 GeV [7], assuming B(t' W{q = 
d, s, b}) = 100% and B(b' -> Wt) = 100% respectively. 
This suggests that fourth generation fermions must in- 
deed be heavy, in support of the compositeness scenario. 
These searches typically have been interpreted under the 
assumptions of m t > — my < M w and negligible mixing of 
the (t', b') states with the two lightest quark generations. 
Such conditions are generally required for the SM4 with 
one Higgs doublet, to account for EW precision data [12]. 

Moreover, when a fourth generation of fermions is em- 
bedded in theories beyond the SM, the large splitting case 
(m t >—my > Mw) and the inverted scenario (rrtf < my) 
have not been excluded. An example of this was given 
recently [5], showing that precision EW data can acco- 
modate m t ' — my > Mw if there arc two Higgs doublets. 
In fact, the compositeness picture emerging from the ad- 
dition of new heavy fermionic degrees of freedom is more 
naturally described at low energies by multi-Higgs theo- 
ries [3-5]. 

This subtle theoretical landscape suggests that there is 
no uniquely interesting set of assumptions under which 
experimental data must be interpreted. From the ex- 
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pcrimental view point, choosing a simple set of assump- 
tions allows straightforward, if narrow, interpretation of 
results. Yet these interpretations may be extended [9] 
to give broader and more general results. This article 
continues the work of our previous Letter [9]; we again 
consider the possible 4th generation flavor-mixing space 
broadly, and apply results of the previous searches by 
CDF to calculate direct limits for arbitrary mixing val- 
ues and for both cases of \mt> — my\ > Mw- This arti- 
cle includes significant unpublished details of the previ- 
ous calculations, an additional recently published third 
dataset [7], and a translation of the mass limits that we 
derive in branching fraction space into mass limits in fla- 
vor mixing space, which allows for direct comparison with 
other experimental flavor mixing constraints. 

The following sections describe the original CDF mea- 
surements, mass limits in branching fraction space, and 
the limits remapped to CKM4 space. 



II. SAMPLES AND STRATEGY 

The CDF Collaboration has published several impor- 
tant limits on fourth-generation quark masses. We ori- 
ent this discussion with a summary of their samples and 
results. Three analyses have been presented: (1) a col- 
lection of events containing at least four jets and a single 
lepton, known as the I + 4j sample[6]; and (2) a col- 
lection of events containing a single lepton, and at least 
five jets (one with a flavor tag), known as the I + 5j 
sample [7]; and (3) a collection of events containing two 
same-charge leptons, two jets (one with a flavor tag) and 
evidence of neutrinos (missing transverse energy) , known 
as the e ± £ ± jb^ T sample [8]. 



A. The £ + 4j sample 

In 4.6 fb _1 of data, CDF searched for t' W{q = 
d, s, b} decays in the mode 

pp -> t't' -t {W -> lu)q(W -> qq')q 
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by requiring a single lepton and at least four jets. The 
data were analyzed by reconstructing the invariant mass 
of the candidate t' and measuring the total energy in the 
event. The event selection used the four jets of high- 
est transverse energy in the event, but did not require 
a flavor signature (fe-tag) on any of the jets, making it 
generally sensitive to t' — > W{q = d,s,b}. Assuming 
B(t' — > W{q = d, s, b}) = 100%, CDF found m' t > 335 
GeV. 

The mass reconstruction used minimum-likelihood fit- 
ting methods that depend upon the particular spectra of 
final state components. If, for example, one half of an 
event decayed as t' — > W(b' — > Wq) giving a WWqWq 
topology (rather than t' — > Wq giving the WqWq topol- 
ogy) , it might satisfy the I + Aj selection criteria, but the 
reconstructed mass distribution for such events would be 
significantly modified by the additional W. Thus, the re- 
sults cannot be trivially applied to topologies other than 
WqWq, even when they are expected contributors to the 
studied sample. We therefore apply the t + Aj results to 
WqWq processes exclusively. 



B. The £ + 5j sample 

In 4.8fb~ 1 of data, CDF searched for b' -> Wt decays 
in the mode 



pp^b'V ^WtWt -> WWbWWb 

-+ (e ± iy)(qq')b(qq')(qq')b 

by requiring at least one lepton, at least five jets (one 
with a b flavor tag) , and missing transverse energy of at 
least 20 GeV. The data were analyzed by examining the 
number of jets in the event and the total scalar transverse 
energy in the event (H T ). Assuming B(b' — > Wt) — 
100%, CDF found m' b > 385 GeV with this signature, 
the highest limit to date. 

As in the i + Aj sample, the analysis uses minimum- 
likelihood fitting methods which requires knowledge of 
the distribution of the signal and background in the Ht 
variable as well as jet multiplicity. Rcinterpretting these 
results to set limits on mixtures of b' — > Wt and t' — > Wb' 
is not possible without knowledge of these distrubutions, 
so we apply the t + 5j results exclusively to the WtWt 
processes. 



The e ± e ± jbfl T Sample 



CDF also searched for b' -> Wt decays in 2.7 fb 1 of 
data, in the same-charge lepton mode 

pp -> b'V -> WtWt -> WWbWWb 

-+ (£ ± ^(qq')b( qq ')(e ± ^b 

by requiring two same-charge leptons, at least two jets (at 
least one with a 6-tag), and missing transverse energy of 
at least 20 GeV. 



Given the small backgrounds, multiple neutrinos and 
large jet multiplicity in the sample, CDF did not recon- 
struct the b' mass, but instead fit the observed jet mul- 
tiplicity to signal and background templates generated 
from simulations. Assuming B(b' -> Wt) = 100%, CDF 
found m' b > 338 GeV. 

Since the l^l^jb^T analysis did not use final-state de- 
pendent fits, these results are process independent — they 
may be applied to any process producing the i^l^jb^T 
signal. For example, t' -> Wb' -> WWt -> WWWb de- 
cays would produce a six- IT, two-6 signature, with higher 
jet multiplicity and larger acceptance to the l^l^jb^T 
sample than the simple four-W, two- 6 signature. In this 
analysis, we therefore apply the l^l^jbT^T results to four- 
W, two-6 processes inclusively. 



D. Analytic Extension of the CDF Results 

The three data samples can be seen as probing a region 
of a two-dimensional interval in branching fraction space 
which, for the classical splitting (t' heavier than &'), is 
specified in the following way. 
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FIG. 1: The flavor-mixing interval overlaid with a table of the 
processes contributing to the corner points 

The topologies of b' and t' decays are determined by 
four branching fractions, two of which are independent: 



B(t' -> Wb') 
B(b' -> Wt) 



B(t' -> W{q : 
B{b' -> W{q 



d,s,b}) 
u,c}) 



The dependence among these quantities, and the pro- 
cesses they represent, is shown in Fig. 1. In this repre- 
sentation, the £ + Aj analysis, l^l^jb^T and £ + 5j anal- 
yses probe the corner (0,1). Each used the assumption 
of an individual contribution from one flavor of fourth- 
generation quark. 

We consider the implications of the CDF results in 
flavor-doublet scenarios, characterized by different as- 
sumed mass splittings and a continuum flavor-mixing 
rates. To extend the interpretations of the published re- 
sults, we use the basic relationship among event yield, 
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cross section and acceptance to reinterpret the observed 
yield limits under different assumptions. This requires 
careful estimation of the relative acceptance between the 
original assumptions, and those of the new interpreta- 
tion. 



III. MASS LIMITS 



B( t'-> Wq ) 



B( t'-» Wq ) 



T 0.8 
«0.6 
0.4 
0.2 



0.2 0.4 0.6 0.8 1 
B( t'-»Wb' ) 




0.2 0.4 0.6 0.8 1 

B( t'^W ) 



FIG. 2: Limits on b' mass from I + 4j data, as a function 
of branching fractions B(t' -> Wb') = 1 - B(t' -> Wq) and 
Bib' — > Wt') = 1 — B(b' — > Wq) for mass structures m t i = 
my + lOOGeV (left) and m t , = my + 50GeV (right). 
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FIG. 3: Limits on b' mass from £ + 5j data, as a function 
of branching fractions B(t' -> Wb') = 1 - B(t' -> VFg) and 
B(6' -> WY) = 1 - B(6' -> Wq). The limits do not depend 
on rn t r — my as t' gives no contribution. 



FIG. 5: Limits on b' mass from e ± £ ± jb$ T , £ + 5j and £ + 4j 
data, as a function of branching fractions B(t' — > Wb') = 
1 - B(t' -> Wg) and B(6' -> Wt') = 1 - B(b' -> Wq) for mass 
structures m t / = m b / + lOOGeV (top) and m t < = my + 50GeV 
(bottom). 
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FIG. 4: Limits on b' mass from t ± ( ± jb^ T data, as a function 
of branching fractions B(t' -> Wb') = 1 - B(t' -> Wq) and 
S(6' — > Wt') = 1 — B(6' — > Wq) for mass structures m t i = 
my + lOOGeV (left) and m t , = my + 50GeV (right). 



In general, the event yield divided by the integrated 
luminosity, N/L, is equal to the cross section times the 
acceptance rate. For a particular process, such as indi- 
vidual b' , this simply tells us that the limit on a cross 
section, given an observed yield limit, is given by 



N 



Ob' 



L ■ e v 



(1) 



where ey is the acceptance rate for the observed process 
within the experimental selection constraints. However, 
we can also consider the case with two contributions if we 
know the relative acceptance rates between the processes, 
e rei , and if the two cross sections are dependent: 



N 



ay + e rel ■ 0- t i(<Jy). 



(2) 
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Here, the dependence of the t' cross section on the 
b' cross section is given, from the next-to-leading order 
cross-section calculations for massive quarks[14] , by the 
mass splitting. Need to add comment on sigma rell 

To probe the full two-dimensional branching fraction 
interval, we must calculate the dependence of the event 
yield on the branching fractions explicitly. As the branch- 
ing fractions to the reconstructed states vary, the accep- 
tances of the processes of interest vary accordingly. Re- 
construction efficiencies can differ as well. Considering 
these effects, we calculate the acceptances of b' and t' in 
the £ + 4j, £ + 5j and £ ± £ ± jb^ T samples. 

The expected signal yield for one process relative to 
another is proportional to relative production rates and 
final-state reconstruction efficiencies. The relative signal 
production rates have two factors: the relative cross sec- 
tion for initial state production and the branching ratios 
of the involved processes (second-order when both sides 
of the event are considered, as here). When several pro- 
cesses contribute to a signal there are multiple terms of 
this form. Fixing the event yield at its observed value, 
we may isolate the cross section we wish to limit, and 
express it in terms of the previously measured limit and 
an effective relative acceptance. The effective relative 
acceptance A includes acceptance terms for all processes 
considered, each scaled by relative cross-section. 

In the £ + Aj case, there are no relative reconstruction 
efficiencies to consider, and we have a simple expression 
for the relative expected yields as a function of (3 b > = 
B(b' -> Wt) and ft, = B(t' -> WV): 

MPb>,Pt>) = (1 - A,') 2 + *rel(l - M 2 (3) 

We use this expression to produce limits on the mass of 
the b' as a function of fourth-generation branching frac- 
tions. The results, shown in the interval introduced be- 
fore, are presented in Fig. 2. 

The £ + 5j case is similar, as we are limited in our 
reinterpretation of the published results, and can only 
write 

Atf v ,fc) = Pl (4) 
with the limits shown in Fig. 3. 

The H^t^jbfir ca se is somewhat more complicated. In 
addition to b' , there are two significant t' processes that 
produce the signal selected for the sample analyzed: 

tf^Wb' -> WWt -> WWWb 
t'^Wb' -> WWc 

Moreover, because of jet multiplicity and the flavor-tag 
requirement, it is necessary to consider relative recon- 
struction efficiencies for each contribution. These factors 
are denoted €nw, where N is the number of interme- 
diate IF-bosons in the process described. Presented in 
Table I, they are estimated from simulated data and are 



TABLE I: Reconstruction efficiencies for t' decays mediated 
by W bosons, relative to the b' reconstruction efficiency of the 
original CDF analysis. 
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calculated relative to the four-VF case considered in the 
original analysis. 

A(p b ,,!3 t ,) = ft,+fi^[(l-/3 b ,fe 4W e cc 

+ 2/3 b /(l - f3 b ,)e 5W e cb + /3 b ,e 6 w£bb] (5) 

The different terms in the t' contribution have different 
reconstruction efficiencies due to jet-flavor tagging, ex- 
pressed by e/i/a where /1/2 is the flavor combination of 
the quarks in the final state of the quark-level t' decay. 
These efficiencies are given by statistics from the raw ef- 
ficiencies of the jet-flavor tag to select the beauty (60%) 
or charm (15%) flavor in a jet: 

£bfc = 1 - (1 - £b) 2 

e c b = e c (l - e b ) + e b (l - e c ) + e b e c 

e cc = l-(f-e c ) 2 

Again, we produce limits on the mass of the b' as a func- 
tion of fourth-generation branching fractions. The results 
are presented in Fig. 4. 

As expected, inspection of these results reveals the 
complementary sensitivities of the CDF analyses. How- 
ever, they are not orthogonal and thus cannot be statis- 
tically combined. In view of this, we produce the best 
limits over the branching fraction interval by choosing 
the stronger of two limits at each point. The combined 
results are presented in Fig. 5. 

IV. CKM4 MAPPING 

To compare results of many different experiments, it is 
useful to decouple the various physical parameters wher- 
ever possible and express the results in terms of the un- 
derlying physical quantities at the theoretical level. In 
this case, these are the elements of the flavor mixing 
matrix, on which there are important experimental con- 
straints. 

The description of SM quark flavor mixing falls to 
the CKM matrix, which, in the presence of a fourth 
generation, is extended to CKM4. It has 16 elements 
parametrized by 4 mixing angles and three irreducible 
CP- violating phases[10]. Present experimental data has 
been brought to bear on the extent to which this could 
be consistent with observation (cf. [10]). 

The ideal case would be to form a clear basis of com- 
parison between theory and experiment, and some of the 
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challenges of doing that are topics of discussion here. The 
primary difficulty is in decoupling the many dependent 
parameters involved, so one aims to formulate analysis 
to this end. Here, we eliminate many degrees of freedom 
by treating light quarks as indistinguishable. The scale 
of light quark mass is relative to the mass difference with 
W, as required for weak decay. This treatment is also 
natural from the perspective of experimental design. 

In the limit of massless particles, the branching frac- 
tions of the previous sections translate directly to the 
squares of the magnitudes of the corresponding CKM4 
elements: 




FIG. 6: Mapping from branching fraction to CKM4 space for 
mass structure m t i = my + 100 GeV. Left, B(t' — > Wb') as 
a function of CKM4 parameters Vaa- Right, £>(&' — > Wt) as a 
function of V43 and V4.4. 
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In the case of massive particles, the partial widths 
carry on this proportionality, but there are also phase- 
space constraints to consider. The branching fraction of 
the massive particle is of course constructed from the 
partial widths of the various possible modes. Thus the 
branching fraction has dependence on the masses of the 
parent and daughter particles. 

This presents a challenge: can we discuss fourth gener- 
ation mass limits in the context of a CKM parametriza- 
tion? Ignoring for the moment the issue of mass depen- 
dance, we can construct a transformation between the 
branching fractions and the CKM elements. 

The plane we have chosen, with axes Bit! — > Wb') = 
1 - B(t' -> Wq) [q = d, s, b] and B(b' -> Wt) = 
1 — B(b' — > Wq) [q = u, c] includes processes with tree- 
level diagrams having vertex factors of V44 and V34 re- 
spectively. Of course, a branching fraction is the quotient 

-^partial/^totah 

At tree-level, the partial widths of the final states are 
the product of phase-space and weak-coupling factors, 
and the absolute square of the tree-level CKM vertex 
factor. 

For the total widths of the denominators, we would like 
to eliminate dependence on the mixing angles between 
the fourth generation and the lighter quarks to minimize 
the number of variables involved in the transformation. 
At tree-level, this seems again possible: 

1. The phase-space factors of matrix elements to the 
various flavors are nearly identical for all light 
quarks, including b when the mother quark is at'. 

2. By weak universality we can eliminate the CKM 
parameters that describe the mixing with lighter 
quarks. 



A. B(t' -> Wb') as a function of V i4 . 

For the numerator: 



r(t',b') 

II^IMI J dp\...)f 



The integral over four-space is meant to signify the full 
matrix clement for the process indicated by superscript, 
excluding only the CKM factor. The matrix element is 
calculated using the software program BRI [11], in a con- 
figuration that models no quark mixing (to decouple the 
CKM suppression). 

For the denominator, we have such a term for each 
generation: 
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In the case of the branching fraction for unmixed weak 
fourth-generation decay, we obtain a simple transforma- 
tion between the branching fraction and the correspond- 
ing CKM vertex factor. 



B. B(b' -> Wt) as a function of V43 and V44. 

For the case when the fourth generation b' mixes with 
the third generation in decay, the situation is slightly 
more complicated. Here we cannot describe the branch- 
ing fraction using only the CKM vertex factor V34, be- 
cause the presence of the kinematically inaccessible (as- 
suming classical splitting) t' nevertheless affects our ex- 
pression for the total width when we apply the weak uni- 
versality constraint to remove the dependence on mixing 
with the lighter generations. 



FIG. 7: Mapping from branching fraction to CKM4 space for 
mass structure m t i = ray + 50 GeV. Left, B(t' — > Wb') as 
a function of CKM4 parameters V43 and V44. Right, B(b' — >■ 
Wt). 




FIG. 8: Limits on b' mass from £ ± £ ± jb$ T , £ + 5j and I + 4j 
data, as a function of CKM4 parameters V43 and V44 for mass 
structures m t i = my +100 GeV (top) and m t i = my +50 GeV 
(bottom). 
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In this simplified model, the transformation has only 
two degrees of freedom associated with the phase space 
integrals. These may be expressed by the ratio of the 
two factors in each equation. This is expected since we 




FIG. 9: Limits on b' mass from t^e^jbfh, £ + 5j and £ + 4j 
data, as a function of CKM4 parameters V43 and V44 for mass 
structures m t i = ray + 100 GeV (top) and m t i = ray + 
50 GeV (bottom) using alternate mappings from branching 
ratio space with varied phase-space assumptions. Compare 
to Fig. 8. 



have reduced the structure of the quark masses to two 
splittings, that between t' and 6', and between b' and t. 

In Fig. 6, we show the branching ratio as a function of 
the two relevant CKM4 parameters, V43 and V44 for the 
mass structure m t i = m b > + 100 GeV. Figure 7 shows 
the same for % = my + 50 GeV. Although strictly 
speaking the W is virtual in this case, its behavior in the 
decay is essentially indistinguishable from the on-shell 
case. Finally, we present limits on the fourth generation 
quarks in the CKM4 space, see Figure 8. 

The transformation is sensitive to the phase space as- 
sumptions, as the amount of available phase space is 
mass-dependent. Thus it is important to consider the 
consistency of this assumption with the limits set. To 
demonstrate this, we perform the transformation with 
phase-space assumptions that span the range of limits 
set. The degree to which the limit contours in CKM 
space change is a measure of the robustness of the trans- 
formation. There is very little change as the mass of the 
b' increases to the the top of the range and beyond. As 
the b' mass is reduced to 300 GeV and below, some appre- 
ciable variation occurs. This variation is shown in Fig. 9 
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be robust over the interval of limits set by these data. 

The region favored by other experimental data is V43 < 
0.3 and V44 > 0.94. Figure 10 shows the limits in this 
region of interest. 



V. CONCLUSIONS 

We reiterate the conclusions of our earlier Letter [9]: 
we find that the CDF data imply limits on my and m t > 
of 290 GeV and greater over the full range of mixing 
scenarios, for two characteristic choices of the t'- b' mass 
splitting: m t ' > my and my > m t i . The inclusion of 
a t' strengthens the previously obtained b' mass limit 
from the i^l^jbJSx sample; in the m t / > my case, by 
up to 10% when < m t ' — my < Mw- Because of its 
range of sensitivity, the addition of the £ + 5j sample, 
while increasing the mass limits in some neighborhoods 
of parameter space, does not significantly improve the 
overall results of this analysis. 

The transformation of the mass limits to CKM space 
is found to be fairly robust despite the variation in the 
phase-space description. Based on electroweak CKM 
limits[10], we can identify a region of branching-fraction 
parameter space that is most interesting, but this does 
not further constrain the range of mass limits set by the 
available data. 



FIG. 10: Limits on b' mass from f ± £ ± jbft T , I + 5j and £ + 4j 
data, as a function of CKM4 parameters V43 and V44 in the 
region of interest, for mass structures m t i = my + 100 GeV 
(top) and m t i = my + 50 GeV (bottom). 



. Inspection of the figures shows the transformation to 



VI. ACKNOWLEDGEMENTS 

We thank Tim Tait, Shaouly Bar-Shalom and Alexan- 
der Lenz for useful conversations. The authors are sup- 
ported by grants from the Department of Energy Office 
of Science and by the Alfred P. Sloan Foundation. 



[1] P.H. Frampton, P.Q. Hung and M. Sher, Phys. Rept. 
330, 263 (2000); B. Holdom, W.S. Hou, T. Hurth, M.L. 
Mangano, S. Sultansoy, G. Unel, presented at Beyond the 
3rd SM generation at the LHC era workshop, Geneva, 
Switzerland, Sep 2008, arXiv:0904.4698 [hep-ph], pub- 
lished in PMC Phys. A3, 4 (2009). 

[2] B. Holdom, Phys. Rev. Lett. 57, 2496 (1986) [Erratum- 
ibid. 58, 177 (1987); W.A. Bardeen, C.T. Hill and M. 
Lindner, Phys. Rev. D41, 1647 (1990); C. Hill, M. Luty 
and E.A. Paschos, Phys. Rev. D43, 3011 (1991); P.Q. 
Hung and G. Isidori Phys. Lett. B402, 122 (1997). 

[3] P.Q. Hung and Chi Xiong, arXiv:091 1.3890 [hep-ph]; 
ibid. arXiv:091 1.3892 [hep-ph]; H. J. He, C. T. Hill 
and T. M. P. Tait, Phys. Rev. D 65, 055006 (2002) 
[arXiv:hep-ph/0108041]. 

[4] M. Hashimoto, V.A. Miransky, arXiv:0912.4453 [hep-ph]. 

[5] M. Hashimoto, arXiv: 100 1.4335 [hep-ph]. 

[6] A. Lister (for the CDF collaboration), pre- 
sented at ICHEP 2008, arXiv:0810.3349 [hep-ex]; 
J. Conway et al, CDF public conference note 
CDF/PUB/TOP/PUBLIC/10110. 



[7] L. Scodellaro (for the CDF collaboration) presented at 
ICHEP 2010. D. Whiteson et al, CDF public conference 
note CDF/PUB/TOP/PUBLIC/10243 
[8] T. Aaltonen et al. (by the CDF Collaboration), 

arXiv:0912.1057 [hep-ex]. 
[9] C.J. Flacco, D. Whiteson, T.M.P. Tait, S. Bar-Shalom, 
Phys. Rev. Lett. 105, 111801 (2010) 
[10] M. Bobrowski, A. Lenz, J. Riedl, J. Rohrwild, Phys. Rev. 

D79, 113006 (2009); 
[11] P. Meade, M. Reese, arXiv:hep-ph/0703031v2 [hep/ph] 
[12] G.D. Kribs, T. Plehn, M. Spannowsky, T.M.P. Tait, 
Phys. Rev. D76, 075016 (2007); ibid. Nucl. Phys. Proc. 
Suppl. 177-178, 241-245 (2008); M.S. Chanowitz, Phys. 
Rev. D79, 113008 (2009); V.A. Novikov, A.N. Rozanov 
and M.I. Vysotsky, arXiv:0904.4570 [hep-ph] and refer- 
ences therein; J. Erler and P. Langacker, arXiv:1003.3211 
[hep-ph]. 

[13] A. Arhib and W.S. Hou, JHEP 0607, 009 (2006); M. 
Bobrowski, A. Lenz, J. Riedl and J. Rohrwild, Phys. 
Rev. D79, 113006 (2009); A. Soni, A.K. Alok, A. Giri, 
R. Mohanta, S. Nandi, arXiv:0807.1971 [hep-ph]; G. 



8 



Eilam, B. Melic and J. Trampetic, Phys. Rev. D80, 
116003 (2009); A. Soni, A.K. Alok, A. Giri, R. Mohanta, 
S. Nandi, arXiv: 1002.0595 [hep-ph]; A.J. Buras et al., 
arXiv:1002.2126 [hep-ph]. 

[14] R. Bonciani, S. Catani, M. L. Mangano, and P. Nason, 
Nucl. Phys. B529, 424 (1998). M. Cacciari, S. Frixione, 
M. L. Mangano, P. Nason, and G. Ridolfi, J. High Energy 
Phys.04 (2004) 068. 

[15] H. J. He, N. Polonsky and S. f. Su, Phys. Rev. D 64, 



053004 (2001) [arXiv:hep-ph/0102144]. 
[16] R. Bonciani, S. Catani, M. L. Mangano, and P. Nason, 

Nucl. Phys. B529, 424 (1998); M. Cacciari, S. Frixione, 

M. L. Mangano, P. Nason, and G. Ridolfi, J. High Energy 

Phys. 04 (2004) 068. 
[17] J. Alwall et al., JHEP 0709, 028 (2007) [arXiv:0706.2334 

[hep-ph]]. 

[18] P. Meade and M. Reece, arXiv:hep-ph/0703031. 



